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Abstract 

We derive the Boltzmann equation for scalar fields using the Schwinger-Keldysh for- 
malism. The focus lies on the derivation of the collision term. We show that the rele- 
vant self-energy diagrams have a factorization property. The collision term assumes the 
Boltzmann-like form of scattering probability times statistical factors for those self-energy 
diagrams which correspond to tree level scattering processes. Our proof covers scattering 
processes with any number of external particles, which come from self-energy diagrams with 
any number of loops. 
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I. INTRODUCTION 



Kinetic theory has proven to be a very successful tool for the description of systems 
out of equilibrium. Important applications can be found in many fields of physics, 
current examples are early cosmology, or the theory of heavy ion collisions, where 
one aims to understand how the matter produced in the collision evolves, and in 
particular whether it thermalizes or not. A very famous kinetic equation is the Boltz- 
mann equation for the evolution of the particle distribution functions in a sufficiently 
dilute system. A first derivation of non-relativistic kinetic theory from field theory 
was given by Kadanoff and Baym Using the Closed-Time-Path formalism, also 
called Schwinger-Keldysh formalism, this derivation becomes simpler and can also be 
extended to relativistic theories, see for example In early treatments 

the collision term of the Boltzmann equation was simply modeled heuristically: con- 
sider each reaction the particle can undergo, compute the probability for this reaction 
by using the free particle cross section, and multiply with the appropriate statistical 
factors, that is the Bose-enhanced or Fermi-suppressed phase space distribution func- 
tions, respectively. The collision term for a scalar particle that can undergo 2-to-2 
scattering, for example, is 

= 4^ / i2n)L, {2t)L, {2n)Lj ^''^'^'^^' + ^' ' ^ ' ^) 

X V^up'k' X [U'fk'{l + fp){l + fk) - (1 + U'){1 + fk')Ufk\. 

We refer to a collision term of this form as a standard Boltzmann collision term. In 
references 

mm 

the collision term of the resulting Boltzmann equation was de- 
rived for a number of specific theories from the right hand side of the Kadanoff-Baym 
equation considering self-energy diagrams with up to two loops. To our knowledge 
the first paper with a computation that went beyond the 2-loop self-energy is 0]. The 
authors explicitly computed a selection of self-energy diagrams with up to three or 
four loops in real, scalar 0^- and 0^-theory, respectively, and managed to bring them 
into a form like (Q). This computation was quite involved and required the assistance 
of a computer system to handle the big number of terms appearing in intermediate 
steps. 

In the present paper we consider a scalar theory with an unspecified non-derivative 
self-interaction. We first show that the self-energies appearing in the collision term 
have a useful factorization property. In the strict on-shell limit a certain class of self- 
energy diagrams then indeed leads to a collision term of the form (0), namely those 
that correspond to tree level scattering diagrams with any number of external parti- 
cles. Not all contributions to the self-energy fit into this picture, which can already 
be expected from the problems that arise when extending the vacuum Cutkosky rules 
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to finite temperature 




In section IH] we give the basics of nonequilibrium field theory as far as required for 
our purpose. We also show briefly how the flow term of the Boltzmann equation is 
obtained from the Schwinger-Dyson equation. This procedure is standard and has 
extensively been covered in literature already, so only the basic steps are given here. 
For more details see for example and references therein. Section IIIII is the main 
part of this paper, where we study the collision term for a real scalar field. A remark 
on the extension to charged scalar fields can be found in section IIVI and in section |3 
we finally discuss our results. 



II. BASICS 

An appropriate framework to study the time evolution of physical quantities in 
nonequilibrium situations is given by the Schwinger-Keldysh or Closed- Time-Path 
formalism 10. 

HQ. 

The basic technical point is that the time variable of all ob- 
jects is defined on a path C on the real axis that leads from — oo to +00, and then 
back to —00. The definition of the Green function is still 

A{x,y) = -i{Tcct>{x)<p\y)), (2) 

but time ordering is performed along the path C here. If we split the path into a "+"- 
branch from —00 to +00 and a "—"-branch going back to —00, we can distinguish 
four real time Green functions, differing by the branches on which the time arguments 
are placed: 

zA++(x, y) = ^A*(x, y) = {T<P{x)4>Ky)) , (3) 
iA+'{x,y) = iA<{x,y) = {(t)^{y)(j){x)) , 
iA~^{x,y) = iA>{x,y) = {(f){x)(l)\y)) , 
iA-{x,y) = iA\x,y) = {f(t){x)(t)\y)). 

The T in the last line denotes anti time ordering. These functions are not independent 
of each other, but are connected via the relation 

A* + A* = A< + A>. (4) 

In addition, one defines the retarded and advanced Green functions: 

A«'^(x,i/) = A*(x,i/)-A<'>(x,y). (5) 

In contrast to ordinary perturbation theory in vacuum, here each internal vertex can 
be either of type "-|-"or of type "— ". The lines between vertices represent the Green 
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functions defined in Q, depending on the types of the incident vertices. For each 
vertex of type "—"there is an additional factor —1. 

The Schwinger- Dyson equation for the contour Green function (j2I) leads to the fol- 
lowing equation of motion: 

(5^ + m') A{x, y) = -S'cix - y) - [ d'z n(x, A(^, y). (6) 

Jc 

The self-energy 11 is defined as i times the sum of all one-particle irreducible two-point 
functions. We switch to index notation and obtain the equation of motion for A^''*: 

{dl + m^) A<'>{x, y) = j d^z (n<'>(x, z)A^{z, y) - n^(x, z)A<^>{z, y)) . (7) 

For simplicity we omit a potential tadpole contribution to the self-energy 11. We 
assume that its only effect is a shift in the mass of the particles which possibly 
introduces a force on the left hand side of the resulting Boltzmann equation, but that 
doesn't change the argumentation concerning the collision term. 
Our goal is an equation that describes the evolution of the phase space density of 
particles, which takes place on a macroscopical scale. In order to separate this from 
the quantum evolution on a microscopical scale, we perform a Wigner transformation. 
For any two-point function we introduce the average coordinate 

X = \{x + y) (8) 
and carry out the Fourier transformation with respect to the relative coordinate: 

A(X,p) = j d\x-y)e'P-^^-y^A{x,y). (9) 

Note that the functions zA'^'^(X, p), called Wigner functions, are real, while zA*(X, p) 
is the complex conjugate of zA*(X, p). The respective self-energies have the same 
properties. 

We apply the Wigner transform to the equation of motion ((Tj). The fact that we are 
dealing with two-point functions is reflected in the appearance of an inflnite series of 
derivatives 

o{(l)}{(2)} = \ ■ 9f - 5« ■ 9(2)) {(1)}{(2)} (10) 
with respect to X and p: 

-p^ + -ip-d+-dA A<^>{X,p) (11) 
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{n<'>(x,p)}{A^(x,p)} - {n^(x,p)}{A<'>(x,p)}). 



Now some approximations are necessary. We assume that the functions of interest 
have a smooth macroscopic behavior, more precisely, we assume that the scale on 
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which the functions change with respect to the average coordinate X is much bigger 
than the microscopical scale set by the de Broglie wavelength of the particles. Con- 
sequentially, the mixed derivative d ■ dp is a. small quantity and kept only up to linear 
order. We furthermore assume that the macroscopical scale is also large compared 
to the particles' Compton wavelength, so that the second order derivative on the left 
hand side of (|TT| is negligible as well. Since the Wigner functions zA^'-^ are real, we 
can extract the real and imaginary part of the equation of motion and find 

{p^ - m^) iA<'> = (n<'>ReA^ + Ren^A<'>) - ^o(n<A> -n>A<),(12) 

{-p ■ d) tA<'> = ^ (n> A< - n<A>) - ^ o (n<'> Re A^ + Re A<'>) . (13) 

The real part has the form of a constraint, the imaginary part has the form of a 
kinetic equation. 

For a free field we can write down the solutions of these equations in the form 

tA<{X,p) = 27r6{p^-m^)sgnipo)no{X,p), (14) 
iA>{X,p) = 27r5(p2-m2)sgn(po)(l + no(X,p)). 

The solutions for the chronological and antichronological Green functions in the free 
case are 

tAl{X,p) = -— + 27r5(p2 _ m2)sgn(po)no(X,p), (15) 

p^ — + 2sgn(po)e 

tAl{X,p) = +27r6{p^ -m^)sgn{po){l + no{X,p)). 

In thermal equilibrium the KMS relation determines uq to be the Bose-Einstein dis- 
tribution, but in a nonequilibrium situation this function is not known a priori. 
The next simplification is a small coupling expansion. Later this will be used for a 
detailed analysis of the collision term, here we need it to get rid of those terms on 
the right hand side which are suppressed by both the mixed derivative d ■ dp and the 
coupling constant. At this point we also demand that our system can be described 
in terms of (quasi-)particles. To this end we assume that the right hand side of the 
constraint equation (fT^ vanishes, turning this equation into a mass-shell condition. 
In the end, equations (fT^ and (fT^ simplify to 

{p^-m^)iA<^> = 0, (16) 

(-p-9)iA<'> = i(n>A< -n<A>), (17) 

and we can make an on-shell ansatz for the Wigner functions: 

iA<{X,p) = 27i6{p^ - m^)[9{po)f+{X,p) + 9{-po){l + f-{X,-p)) 

iA>{X,p) = 2n5ip' - m') \9{po){l + f+{X,p)) + e{-po)f-{X, -p) 



5 



Spectral sum rules that follow from the basic commutator relations for the scalar field 
operators and make a connection between zA< and iA-* ensure that this ansatz is 
consistent. By comparison with the equilibrium functions and also by inserting the 
ansatz (fTHj) into the expressions for the expectation values of current 

nX) = 2 I ^K^A<(X,p) (19) 

and energy-momentum 

T'-'iX) = J ^pV^A<(X,p), (20) 

we finally can identify /+ and /_ with the phase space densities of particles and anti- 
particles, respectively. In the case of real scalar fields there is an additional relation, 
iA'^i^—p) = iA^{p), which leads to /+ = /_.. 

The constraint equation is satisfied identically with this ansatz, and all that is left over 
is the kinetic equation (fTTj) . We insert (fTH|) and integrate over positive frequencies. 
The resulting equation reproduces the flow term of a relativistic Boltzmann equation 
for particles with phase space density /+: 

{d, + ^- UiX,p) = C{X,p). (21) 

Integration over negative momenta results in a similar equation for the corresponding 
antiparticles. The most important part for us is the right hand side: the collision 
term so far is 

C{X,p) = ^(iU>{X,p)zA<{X,p)-zU<{X,p)zA>{X,p)y (22) 

In the next section we try to re-express this in terms of particle densities and scattering 
amplitudes. 

III. COLLISION TERM 

Now we come to the main part of this paper, the collision term. We make a pertur- 
bative expansion of the self-energies in C and try to bring it to a form resembling the 
collision term of a Boltzmann equation like shown in the introduction. The collision 
term is local in our approximation, so we drop the argument X from now on in order 
to simplify the notation. 

A. Self-energy 

In a perturbative expansion the self-energy in^(p) = i]l~^~{p) is expressed as the 
sum of all amputated one-particle irreducible graphs with momentum p entering at a 
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FIG. 1: A graph with only dh'ect connec- 
tions between "+"- and "—"-base (type 1). 
There is a unique complete cut (dashed 
line). 



FIG. 2: A graph with indirect connec- 
tions between "+"- and "—"-base (type 2). 
Two different complete cuts are possible 
(dashed lines). 



"+" -vertex and leaving at a "—"-vertex. We first classify the graphs in this expansion 
in the following way: take any graph and imagine all lines connecting a "+" -vertex 
with a "—"-vertex were cut. This would split the graph into a number of connected 
subgraphs, each containing only "+"- or "—"-vertices, respectively. We call these 
subgraphs clusters. Obviously there are at least two clusters, namely one which is 
connected to the incoming line, called "+"-base, and one connected to the outgoing 
line, called "—"-base. Now we can distinguish two types of graphs: 

1. Graphs of type 1 have only direct connections between the "+"- and the "— "- 
base, i.e. there are no paths leading from the "+"-base to the "—"-base via 
some other clusters (see Fig. Q). The simplest example for this are graphs with 
only two clusters. 

2. Graphs of type 2 have connections between the "+"-base and the "—"-base via 
other clusters (see Fig. E}. 

There is another way to classify the graphs contributing to ill~^~{p): try to divide a 
given graph into two parts, one attached to the incoming line, the other part attached 
to the outgoing line, by cutting a suitable set of "H — "-lines. Only cut lines which 
are attached to a "+"-vertex in the first part and to a "—"-vertex in the second part, 
and which are effective in separating the two parts. This works with every graph, 
since the incoming line ends in a "+" -vertex, while the outgoing line starts from a 
"—"-vertex. We call this a complete cut. It is not hard to see that a graph belongs 
to type 1 if and only if there is exactly one such complete cut. In this case the cut 
lines are exactly the ones connecting the "+"-base with the "—"-base. Otherwise the 
graph belongs to type 2. 
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Unique complete cut 



Let us first concentrate on the graphs of type 1, i.e. graphs which have a unique 
complete cut. We further classify these diagrams according to the number of cut 
lines: 

oo 

^nJ,-ei)(p) = E^n+-(p). (23) 

n=2 

Now we claim that 

- w = ^ / S ■ ■ ■ Sf'^'^^'''^" (2*) 

X -^1, ■ ■ ■ , -kn) iK'ih) ■ ■ ■ i^o'ikn) A^^i)(-P, ki,...,kn). 

Symbolically we can write this equation as 




The infinite sum of diagrams that contribute to the self-energy on the left hand 
side factorizes into two other infinite sums which are identified as the diagrammatic 
expansions of certain n-point functions. The totally symmetric function is 
the scattering amplitude belonging to the amputated, connected, out-of-equilibrium 
(n + l)-point function; positive momenta are entering, and all external momenta are 
attached to "+" -vertices. Basically it is the Fourier transform 




of the out-of-equilibrium {n + l)-point function (we use the same symbol for the 
function and its Fourier transform) with time ordered fields 

A+)(a;i,...,a;„+i) = {T(j){xi) . . . <j){xn+i))^^^^, (26) 

however with the additional restrictions that the graphs contributing to this func- 
tion must not have any corrections on the p-line, and for each ki there must be a 
path connecting it to p which only contains "+" -vertices. The derivative operators 
□ = (9^ + remove the external legs and we explicitly take out the momentum 
conservation 5-function, 

A^i)(?i, . . . , Qn+i) = (27r)^(5^(?i + . . . + qn+i)iM^^^{qu . . . , g„+i), (27) 

so that has the form of a scattering amplitude. The functions and A^^ 

are defined analogously with anti-time ordering. 
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In order to prove (j24p . we have to show two things: first, the set of diagrams appearing 
in the perturbative expansion on the left hand side is the same as the set of diagrams 
one finds on the right hand side. Here we don't care about statistical factors of 
diagrams or how often they appear. This is done in a second step: it may be possible 
to obtain the same diagram by putting together several different contributions from 
and A~. The number of such combinations together with the factor l/n\ and 
the symmetry factors of these contributions must match the symmetry factor of the 
resulting diagram on the left hand side. 

Diagrams 

Any graph G contributing to 11^" has a unique complete cut consisting of n lines. 
The free propagators ^Aq"" associated to the cut lines and the integrals over their 
momenta, named fci, . . . , are taken out of G and written down explicitly on the 
right hand side of (j2H)- The remainders of G to the "left" and to the "right" of the 
cut are called G^ and G^, respectively. Obviously G^ is a Feynman diagram with 
n + 1 external momenta attached to "+" -vertices. We took away the fcj-lines, so it 
is amputated. Since G is one-particle irreducible, there can be no corrections on the 
p-line, while there may be corrections on the fcj-lines. Because G is of type 1, the cut 
lines are exactly the ones that connect the "+"-base to the "—"-base, so all external 
momenta in G^ are attached to the "-|-"-base and thus are connected with each other 
via paths that only include "+" -vertices. This shows that G^ is a Feynman diagram 
that contributes to A^-^, and likewise G~ is a part of A^-j. 

Conversely, take any contributions G^ form A^-^ and G~ from A^^. Together with 
the free propagators iA^~ and the integrals over their momenta they make up an 
amputated Feynman diagram G with a unique complete cut and with momentum p 
entering at a "+" -vertex and leaving at a "—"-vertex. Since there are at least two 
fcj-lines, the only way for G not to be one-particle irreducible would be to consist of 
a G"*" or a G~ which can be split into two parts, one connected to the p-line and one 
connected to the fcj-lines, by cutting a single line. But this is impossible, because G^ 
and G^ don't have corrections on the p-line. So the diagram G contributes to the 
diagrammatic expansion of n+~, which eventually proves that the same diagrams are 
appearing on both sides of (j^ . 

A comment about the momentum conservation 5-function is in order here. If one 
performs as many integrals over internal momenta as possible in the diagrams G"*" 
and G~, in both cases a (27r)^5^(p — ki — . . . — kn) is left over which is not part of 
the scattering amplitude. One appears explicitly in (j21I), the other one reduces to 
(27r)'^5'^(0) and is dropped according to the usual definition of the self-energy. 
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Numbers 



Let us first have a closer look at the symmetry factor of a diagram G contributing to 
n+~. The symmetry factor of a diagram is the order of the graph's symmetry group, 
which contains all permutations of lines and vertices that do not alter the diagram. 
G has a unique cut, and accordingly we can distinguish lines and vertices in the part 
G"*" left of the cut from lines and vertices in the part G~ to the right of the cut, and 
they all are topologically different from the cut lines. This means that there are no 
symmetries exchanging lines or vertices in G^ with those in G~ , and neither are there 
symmetries that exchange cut lines with uncut ones. Therefore the symmetry group 
of G is the direct product of three groups: the symmetry group of G^, the symmetry 
group of G~ , and the symmetry group S of the cut lines (in the context of G, i.e. 
permutations of cut lines that do not change G). We call the orders of these groups 
s_ and s, respectively, and so the symmetry factor of G is l/(s_|_ss_). In fact, 
things are a bit more complicated, since in graphs with tadpole-like structures, i.e. 
where both ends of some lines are attached to the same vertex, the above definition 
of the symmetry group is too narrow. However, such lines will never be cut, and 
therefore their symmetry properties are part of the symmetry groups of G~^ or G~ 
and don't change this discussion. 

Given a graph G contributing to , many different combinations of a G^ G 
and a G^ G produce this G? First of all, the topologies of G"*" and G~ are com- 
pletely determined by G and its cut. But there are several possibilities to name their 
external lines, which determines how they are put together. Some namings result in 
topologies different from G, so they contribute to a different graph and are not rele- 
vant here. The rest of the namings corresponds exactly to the inequivalent renamings 
of the cut lines fc, of G, where inequivalent means that they are not symmetries of 
G. So the number of combinations of a G~^ and a G^ that produce G is equal to the 
number of inequivalent renamings or permutations of the cut lines in G. 
The group P of all permutations cxj of cut lines in G, i.e. all possibilities of renaming 
them, has the order n! If we build the right cosets of P relative to S, 

Sai,Sa2, . . . ,San\, (28) 

then two sets Sai and Saj are either identical or disjoint 0|. So in fact there are 
only / different sets 

Sai^,...,Sai_^, (29) 

each of which contains s elements. Since all elements within one Sai belong to 
equivalent permutations of cut lines, the number of inequivalent permutations must 
be equal to /. But the step from ipHjl to only removes redundant elements and 
leaves their total number unchanged: n\ = f ■ s. Together with the symmetry factors 
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FIG. 3: Two 2-loop contributions to the self-energy of real scalar c/i^-theory with a unique 
complete cut comprising n = 3 lines. The symmetry factors of these diagrams are 1/2 and 
1, respectively. 



of the parts and G the complete combinatorial factor on the right hand side 
then is 

4^ ^ ^ (30) 

n! s+S- 

and exactly matches the symmetry factor of G. 



Example 

As an example we consider two 2-loop contributions to the self-energy of a real scalar 
theory with a 0^-interaction, displayed in Fig. |H1 There is only one possible topology 
for the parts G^ and G", namely a 2-to-2-scattering diagram. There are three possi- 
bilities for naming the external lines of this diagram, so that the relevant contribution 
is 



fcl 


1 — fca 1 


ki 


fc2 + 


— fci + 




hi 


— hi 


1 k2 



X 



ki 1 


1 k2 1 


kz 


k2 


+ fci 


+ fci 


hi 1 


kz 1 


k2 



■ (31) 



From the resulting nine graphs, three correspond to the inequivalent ways of naming 
the cut lines in diagram Fig. IHfa), while the remaining six graphs belong to diagram 
Fig.Hb). 



No unique cut 



Next we have to deal with those contributions to the self-energy with more than 
one complete cut (type 2). In this case we decide to use the cut which makes the 
left portion G^ of the diagram as small as possible. This means to cut all lines 
that leave the "+"-base, except when they are not efficient in separating the two 
parts. In the example of Fig. |21 this is the left cut. This prescription leads to an 
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unambiguously defined cut and we can repeat all the arguments already used above. 
The only difference concerns the nature of the right hand portion of the diagrams: 
the contributions to A^^-, must have the property that at least one of the external 
/cj-lines is connected to the outgoing p-line only via paths which contain at least one 
"+" -vertex. Note that and A^(2) don't have any graphs in common. 



B. Collision term 



For real fields the n-point functions A"*" and A are related by complex conjugation: 

A-{h,...,kn) = A+*{-h,...,-kn). (32) 

By referring to the corresponding diagrammatic representation we can check that this 
holds for the restricted functions A^^ and A^^, too. After sending ki —ki the full 
collision term can thus be written 

-'^^t^J§^.-§y. + fa + . . . + fcj (33) 

a=l,2 n=2 -J \ I \ I 



X |;V/l + ^(p,A;i,...,A;„)A<+) (p, fci, . . . , A;„)iA>(A;i) . . . iA>(A;„)iA>(p) 

- M'l^*{-p, -ki, -kn)M^^^{-p, -ki, -kn)iA<{ki) . . . iA<(A;„)iA<(p)|. 

The sum over a = 1,2 represents the contributions of self-energy diagrams of type 1 
and 2, respectively. We emphasize that this factorization is exact. The form (fT^ for 
the Wigner functions suggests that (j^Hj) is an expansion in the particle phase space 
density. But this is not the case, since the distribution function is also contained 
in the scattering matrix. A true expansion in the particle densities has been done 
in which shares some technical similarities with our analysis. 
In order to proceed towards the standard collision term of a Boltzmann equation, 
we still have to overcome three difficulties. First, a matrix element squared can 
obviously only be obtained for type 1 diagrams, where a = 1. Second, the propagators 
corresponding to the cut lines are free ones, while the propagator for the p-line is a full 
one. In the Boltzmann equation, however, all distribution functions are expected to be 
of the same type. And third, if we want to combine the scattering matrices of the last 
two lines, they have to be symmetric with respect to inverting all momenta. The last 
point is satisfied if we assume that our system is symmetric under CP. Note that this 
refers to both the interaction and the initial conditions, since if the evolution starts 
with initial particle distributions / that are not CP-even, then the matrix elements 
won't be CP-even either, even if the interaction conserves CP. The first two points 
force us to restrict ourselves to the tree level, since then the difference between free 
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and full propagators vanishes. In addition, the scattering diagrams include no internal 
"—"-vertices: diagrams with "—"-vertices have several "H — "-lines which cannot all 
be on-shell in a tree diagram, and according to ()14|) these diagrams vanish. 
The fact that we have to exclude quantum corrections is not unexpected. After all, 
the Boltzmann equation is a classical equation, so our derivation definitely has to 
break down at some level when adding quantum effects. As we have seen, this break- 
down takes place immediately beyond the classical level: we can recover a standard 
Boltzmann collision term of the form (P) only by the restriction to classical processes. 
Now we can simplify the collision term to 

(34) 

X |A1+)(p,A;i,...,/c„)|^ 

X (^iA>{ki) . ..iA>{kn)iA>{p) - iA<{ki) . . .'iA<{kn)iA<{p)y 

The final step is to use the on-shell ansatz (fTHj) for the propagators and to perform the 
integrals over the zero-components of the momenta. The ansatz for iA^{ki) contains 
two 5-functions, one that corresponds to positive energies and one that corresponds 
to negative energies. In the case of a negative energy, we additionally invert the 
corresponding spatial momentum. This way we obtain contributions with all possible 
combinations of particles ki to kn either going into or coming out of the scattering. 
Since the scattering amplitude is totally symmetric with respect to the order of its 
arguments, all particles going in are exchangeable, and so are all the particles that 
come out of the scattering. Thus we can group those terms that have the same 
number of ingoing and outgoing particles. For j outgoing particles (besides p) , there 
are 

N{j) = -y^^ = i'") (35) 
identical terms. The complete collision term is 

X (27r)^(5^(p + ki + ... + kj-kj+i- kn) 

^1) • • • 5 ^ji ~kj^i, . . . , —kn) I 

"(l + /l)...(l + /,)/,+l.../„(l + /p) 

-/i.../,(l + /,+i). ..(! + /„)/, 

where = f{ki), uJi = {kf + m^)^/'^, and all four momenta are on-shell: ki = (tUj, ki). 
The collision term consists of two parts, according to the two terms in square brackets. 



X 
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referred to as gain and loss term, respectively. They describe the increase or decrease 
of the density of particles with momentum p in the plasma due to the scattering. 
Typically, several contributions vanish because of kinematical restrictions, for example 
the one with all particles going in. How this works in detail depends on the type of the 
interaction and the particle masses, see also the example in the next section. Since 
in the kinetic equation 1)2111 we only integrate over positive po, terms proportional to 
the negative energy po = —Up were dropped here. 



C. Example: (/)'^-theory at the 2-loop level 



As an example consider a real, massive scalar theory with a X(f)^/3\ self-interaction. 
The relevant 2-loop self-energy diagrams are shown in Fig. IHl The cuts in these 
diagrams comprise three lines, so the corresponding scattering processes in the Boltz- 
mann collision term will have four external lines, and at tree level are shown in Fig. |3] 
The matrix element for these diagrams is 



(p + k2) 



+ 



(p + h) 



+ 



ip+ h 



(37) 



The collision term can then be read off from (jHBj) . In principle there are 4 contribu- 
tions, corresponding to 0, 1, 2 or 3 particles coming out together with p, but due to 
kinematic restrictions only 2-to-2 scattering can occur: 



C{p) 



TT 



■S{po - UJp) 



(27r)32u;i (27r)32a;2 {27r)^2uj3 



(3^ 



X (27r) V(p + k^-k2- ks) \M^ip, ki, -k2, -fcs) r 

;i + /,)(1 + /i)/2/3 - fMl + /2)(1 + fs) 



One could have obtained this result directly from the self-energy computed with the 
OTP Feynman rules. 



2n+-(p) 



d'^q d'^k 
(27r)4 (27r)4 

-iA'+(g)iA++(p + q)iA^'{p + q + k)iA'^{k)iA''{p + q) 



(39) 



+ zA++(p + g)zA+"(p + q + k)iA~~{q + k)iA~^ {q)iA^^ (k) 



but already for these comparatively simple diagrams this is not trivial, and it becomes 
rather involved for self-energies with more loops 0]. 
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k2- 



k-i- 



-ki k2- 




ki ky ^fci 




FIG. 4: The tree diagrams contributing to the scattering matrix for four external particles. 



IV. COMPLEX SCALAR FIELD 

In the case of a complex scalar field, when there are particles and antiparticles, we 
can adopt the general line of argumentation from the real case, but we have to make 
adjustments at a few points. 

First of all, charge conservation constrains the cut: the net charge flow across the cut 
must be equal to the net charge flow through the whole self-energy diagram. Thus 
the cut must comprise 2n + 1 lines, n + 1 leading from the left to the right, and n 
leading back from the right to the left. Then the analog to the basic claim for 
complex fields is 

-m+-in\ = ^ f ^ ^'^"+^ f ^ ^ (AQ) 

" {n + l)\n\ J (27r)4 '■■ (27r)4 J (27r)4 {2n)^ ^ ' 

X (27r)^5^(p + gi + . . . + g„ - fci - . . . - k^+i) 
X A^+^(-A;i, . . . , gi, . . . , g„) 

X zA+-(A;i) . . . zA+-(A;„+i)zAo +(gi) . . .iAo+(g„) 
X -gi, . . . , -g„; /ci, . . . , K+i). 

The proof runs almost exactly like in the real case, the only difference is that we have 
to distinguish the two sets of cut lines. Lines running from the left to the right cannot 
be interchanged with lines running from the right to the left, because they transport 
charge in different directions. This leads to the factor + l)!n! that represents 
one group of n + 1 lines carrying charge to the right and another group of n lines 
carrying charge to the left. 

The final result is similar to ()36p. but both gain and loss term now consist of two parts. 
In one part the distribution functions represent a charge coming out of the reaction, 
where we have all combinations of these functions being either (1 + /+) for a particle 
coming out, or /_ for an antiparticle going in. The other part contains distribution 
functions that represent a charge going into the reaction, that is all combinations 
with either /+ for a particle going in or a (1 + /_) for an antiparticle coming out. 
It is straightforward to generalize to the case of several scalar particle species in a 
similar way. Each species corresponds to a group of lines or distribution functions. 
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respectively, where members of different groups cannot be interchanged with each 
other. 

V. DISCUSSION 

We start with the equation of motion for the out-of-equihbrium Green function for a 
real scalar field. Using gradient expansion, on-shell approximation and small coupling 
expansion the left hand side leads to the flow term of a Boltzmann equation in a well 
known way. We make a perturbative expansion of the self-energies on the right hand 
side and give a diagrammatic proof that the relevant self-energies have a factorization 
property, and that for a certain type of contributions to the self-energy the right 
hand side can be rewritten as the matrix element squared of scattering processes 
times the distribution functions of the particles involved in the reactions. This works 
for those self-energy contributions that correspond to tree level scattering diagrams. 
The scattering processes obtained include any number of external particles and thus 
correspond to self-energy diagrams with any number of loops. 

If effects beyond the classical level are required, it is in general not correct to simply 
compute the scattering probabilities including quantum corrections and plug them 
into a standard Boltzmann collision term. In specific situations it is possible to fit 
quantum corrections into the picture, however at some cost. In [15] for example, a 
model of scalar quarks and gluons was examined and a standard Boltzmann collision 
term was constructed from the self-energy up to two loops. In order to be able to do 
this, the scattering matrices were not computed simply by following the CTP rules, 
but for some diagrams Feynman propagators had to be used for some internal lines. 
Furthermore, kinematical arguments had to be used to get rid of some diagrams that 
didn't match into the picture. The approach presented there required an explicit 
study of each diagram and can not be generalized to arbitrary diagrams or theories. 
The simple structure (J]) cannot be obtained beyond the classical level for two reasons. 
First, there are the self-energy diagrams that allow several complete cuts. This is 
precisely the type of diagrams that causes problems in the extension of the vacuum 
Cutkosky rules to finite temperature |^ . Second, there is the difference between 
the free propagators for the cut fcj-lines and the full propagator for the p-line. Since 
we know that the scattering diagrams we obtain must not contain corrections on the 
p-lines but do have corrections on the fcj-lines, it seems tempting to collect these 
corrections and attribute them to the free propagators, thus promoting them to full 
ones. In effect this means that we would have to do the perturbative analysis of the 
self-energy in terms of resummed and thus full propagators. As a consequence, all 
diagrams contributing to the self-energy must not have corrections to already full 
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lines. Then the available self-energy diagrams are not sufficient to be rewritten as the 
square of a matrix element, however: certain diagrams are missing. For instance, the 
2- loop diagram in Fig. Efb) provides only the interference terms in (jHT|) . the squares 
of the individual amplitudes can only come from Fig. Efa), which has a forbidden 
correction. 

This shows that there is no consistent and systematic way to obtain a standard Boltz- 
mann collision term beyond the tree level. At the tree level, a standard Boltzmann 
collision term is found that includes scattering processes with any number of external 
particles. 
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